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Abstract: We present a new dataset analysis of neutrino experiment at the Baksan Un- 
derground Scintillator Telescope with muon energy threshold about 1 GeV for the longest 
exposure time toward the Sun. In search for a signal from self-annihilations of dark matter 
particles in the center of the Sun we use an update sample of upward through-going muons 
for 24.12 years of live time. The new results are consistent with previous ones [jjj in that 
of no observable excess has been found in measured muons relative to expected background 
from neutrinos of atmospheric origin. We present an improved data analysis procedure and 
describe it in detail. We set the 90% C.L. new upper limits on expected neutrino and muon 
fluxes from dark matter annihilations in the Sun, on the corresponding annihilation rates and 
cross sections of their elastic scattering off proton. 
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1. Introduction 

Evidence for a huge amount of missing (that is dark) mass in the Universe comes from 
observed gravitational effects in astrometry, as it was first discovered by F.Zwicky The 
most impressive hints come from measurements of rotation curves of spiral galaxies and 
gravitational lensing (see Refs. ||, 0], ||, P] for review) including reconstructed view of the 
merging process of bullet clusters 0| , § • Recent high precision measurements of microwave 
background anisotropy [1C] continue to indicate that the total energy density contains only 



about 4.6% of baryonic matter while other parts come from non ordinary forms of gravitating 
matter: dark energy, about 72%, and dark matter (DM), about 23%. The best fit of the 
results of observational astronomy favours the flat AC DM model in cosmology which implies 
existence of the collisionless nonrelativistic DM in an extension of the Standard Model (SM) of 
particle physics. These facts generate numerous multi-wavelength searches for relic particles 
through their possible interactions with ordinary matter by direct and indirect methods. 
Indirect search for dark matter through neutrino channel is accessible within dataset analysis 
of neutrino telescopes in their regular observations of local sources like the Sun where DM 
could be gravitationally trapped and further accumulated for the time of solar system age. 

Up-to now there are no hints on excess of neutrino events in the direction towards the Sun 
as compared with expected background of atmospheric neutrinos at all neutrino telescopes. 
However, these observations allow to set upper limits on properties of dark matter particles, in 
particular, upper limits on their annihilation rates in the Sun which under certain assumptions 
can be translated into upper limits on the cross sections of elastic scattering of dark matter 
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on nucleons. Given the Sun chemical composition with approximately 73% of hydrogen the 
spin-dependent elastic cross section of dark matter particles on proton appears to be one of 
the most sensitive quantities in these searches as it can be seen from the latest results from 
the SuperKamiokande gl|, the IceCube [||] and the ANTARES ||, 0] collaborations. 



The sensitivity of neutrino telescopes strongly depends on their muon energy threshold 
Ej^ 1 , especially in probing for light dark matter particles. The Baksan Underground Scintilla- 
tor Telescope fl5| (further referred to as the Baksan or the BUST) has muon energy threshold 
around 1 GeV allowing for searches for neutrino signal from annihilation of relatively light 
dark matter in the Sun. Corresponding minimal mass which can be probed at the Baksan 
is about 10 GeV. This is intrinsic challenge in light of recent claimed spectral peculiarities 
obtained in low background detectors of nuclear recoil reactions [16|, [18], [|H|, [Q and in 



cosmic ray experiments of 7, e~e + and pp measurements (see Ref. 21] for review). These 
peculiarities have been interpreted by many authors as an evidence on existence of light dark 
matter. 

Here we present updated results of the Baksan experiment with the statistics twice as 
compared to the previous ones [|| and with new analysis improved in several ways. First of 
all, current knowledge of neutrino properties has been fully taken into account. Besides, the 
previous Baksan results Q, |2^| , p3| [ were obtained in the framework of a specific model - 
MSSM (Minimal Supersymmetric Standard Model) with neutralino as a dark matter candi- 
date (see Ref. for a review). Performing a scan over its parameter space but fixing the 
mass of dark matter particle conservative upper limits on muon flux and the rate of dark 
matter annihilation in the Sun have been obtained. In the present work rather then studying 
a specific theory we follow a present day approach which is applicable to a wider class of 
models. Namely, one can make an assumption about dominating annihilation channel for 
dark matter particles in the Sun and then set the limit on the corresponding annihilation 
rate. This approach allows us to compare the Baksan limits with the results of other neutrino 
experiments and at the same time it is to a certain extent model independent: having the 
limits on possible annihilation channels allows one to apply them to a specific model. In our 
analysis we consider annihilations of dark matter particles into bb, t + t~ and W + W~. We 
use our own C code to simulate neutrino propagation from the point of production in the Sun 
to the detector level and compare our results with those obtained with the help of WimpSim 



package |25|, g6[ . 

The paper is organized as follows. In Sec. |2| main parameters of the Baksan Underground 
Neutrino Telescope and levels of neutrino events selection are presented. In Sec. || we describe 
the angular analysis of the upward going muon dataset. In Sec. || we describe in details 
numerical simulations of neutrino propagation and muon signal in the telescope. The results 
and discussion are presented in Sec. || Finally, Sec. |6| contains our conclusions. 

2. Experiment, triggers and sample of upward through- going muons 

The Baksan Underground Scintillator Telescope is a well known aged neutrino telescope 
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located in one of two underground tunnels of the Baksan Neutrino Observatory. The BUST 
is operating cosmic rays detector with 4-pi geometry for arrival penetrating charged particles 
coming from up- and down hemispheres. The duration of continuous measurements covers 
34 years since December of 1978. Main parameters of the telescope have been presented in 
details in Refs. ||15j| . f27fl . Separation of arrival directions between up and down hemispheres 
is made by time-of- flight method with time resolution 5 ns |28[ |. It was shown (2^] that in 
95% of events values of inverse reconstructed particle velocity lie in range of 0.7 -j- 1.3 
for single downward going muons. Such interval but with negative sign of velocity is used to 
select upward going muons generated by neutrino interactions in down hemisphere. 

The telescope is located at the altitude 1700 m above see level in the Baksan valley of 
the North Caucause and at the depth 850 g/cm 2 under the mountain Andyrchi where the 
flux of atmospheric downgoing muons is reduced by factor of about 5000, but it is still higher 
than upgoing muons flux by six orders of magnitude. Trajectories of penetrating particles 
are reconstructed using the positions of hit tanks, which represent together a system of 3,150 
liquid scintillation counters of standard type (70 cm x 70 cm x 30 cm) in configuration of 
parallelepiped (17 mx 17 mx 11 m). The counters entirely cover all its sides and two horizontal 
planes inside at the distances 3.6 m and 7.2 m from the bottom. The planes are separated from 
each other by concrete absorber of radiation length 160 g/cm 2 . The configuration provides 
about 1.5° of muon angular resolution. 

There are two hardware triggers used to select upward going muons [27]. They reduce 
initial rate of downward going muons approximately by factor of 10 3 . Trigger I covers the 
zenith angle range 95° -j- 180° while trigger II selects horizontal muons in the angular range 
80° -7- 100°. The hardware trigger efficiency of 99% has been measured with the flux of 
atmospheric muons [27, 28]. These two triggers give about 1,800 events per day for further 
processing. In the year 2000 the telescope data acquisition system was upgraded and allowed 
to simplify the trigger system down to one general trigger with the rate of 17 Hz. All raw 
information is then undergone further selection using off-line reconstruction code. 

Selection criteria of neutrino events in the Baksan neutrino experiment has not been 
changed since the first data analyses [23], [^?j], [29], [3C[. An off-line program reconstructs 
events and selects those of them which satisfy two simple requirements: presence of only one 
reconstructed track of penetrating particle and negative value of measured velocity f3. In 
addition, all events with negative values of f3 have been scanned by the eyes to check possible 
misinterpretation and there were no observed events within range of —0.5 < /? < 0. Also it is 
required that each trajectory has an enter point lower than exit point on a range not less than 
the tank size. For trajectories crossing only two scintillator planes (sample of trigger II) we 
have excluded tracks having azimuthal angles 0° < < 180° from the direction of minimum 
shallow depth to reduce background from downward going atmospheric muons scattered at 
large angles in the rock. 

The data used for the present analysis have been collected from December of 1978 till 
November of 2009. In total, during 211275 hours of live time (l.t.) which corresponds to 
24.12 years, 1700 events surviving the cuts described above have been collected. 
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Additional cuts have been applied in order to reject upward-going trajectories that could 
be mimicked by downward-going atmospheric muon interactions or multiple muons. i) The 
muon trajectory should have entry and exit points. Thus, stopping muons and neutrino 
interactions inside the detector are excluded, ii) Muon range inside detector must be larger 
than 500 g/cm 2 . These criteria cut off particles with energy below ~ 1 GeV. hi) In the sample 
of events which passed trigger II, those with entry or exit points closer than 1.5 m to the 
plane edge are excluded. This cut reduces background from atmospheric muons interactions, 
iv) Finally, the events with f5 in the range of —1.3 -. — 0.7 only have been accepted. 

In total 1255 events survive all cuts and are used in our further analysis. In Fig. |l] 
(left) we show the rate of collected measured upward through-going muons during all years 
of observation. 

3. The Sun survey during three decades 

In search for neutrinos from the dark matter annihilations in the Sun, we analyse distribution 
of measured events as a function of cos ^-Sun where ^-Sun is the angle between upward 
incoming muons and the Sun position. In Fig. [l] (right) the obtained cosine distribution is 
presented. The mean rate per bin is shown by blue line. We estimate the background expected 
from atmospheric neutrinos directly from real data using shifted (false) Sun positions. Here 
we follow our previous analysis Q, [^2| where it was shown that this method is compatible 
with Monte Carlo (MC) simulations of the detector response. Mean energy of simulated 
atmospheric neutrinos, which produce muons with energy larger 1 GeV to cross the Baksan 
telescope, was about 50 GeV. From comparison of data sample for 21.15 years of l.t. |3(| and 
collected MC statistics larger than real data taking by factor of 22 (i.e., in total 460 years) it 
was found that the ratio of observed total number of events to expected one without neutrino 
oscillations is 0.87 ± 0.03(stat.) ± 0.05(syst.) ± 0.15(theor.). The details have been presented 
in Ref. f(J. 

At the location of the Baksan telescope (43,16°N and 42,41°E) the Sun is seen in average 
about half of time per day during a year in both hemispheres; in Fig. |2| (left) we show solar 
zenith 9 angular distributions reconstructed from data time information (red points with 
error bars) and calculated using the Positional Astronomy Library |3lJ (blue histogram). 
The reconstructed values are obtained from a sample of 1255 selected upward going muons for 
24.12 years of live time. The calculated distribution is normalized to this number. Performing 
the chi-square goodness-of-fit test we obtain a value of x 2 /d.o.f. = 10.38/19. It means our 
multi-years measurements reproduce the Sun full-year passing track with a good accuracy. 
That is the base for our further analysis. 

Potentially, the neutrinos from dark matter annihilations in the Sun can be found at 
night time when the Sun is below horizon. We compare zenith distributions of neutrino 
events coming at nightly (N) and daily (D) time. Corresponding day-to-night asymmetry, 
is shown in Fig. § (left) and no significant difference between these data samples has been 
observed. For "night" neutrinos, we compare cosine distribution in Fig. [| (left) and integrated 
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angular distribution in Fig. ^ (right) with measured background which has been obtained 
from averaged distribution of six cases of false Sun positions shifted along the ecliptic. In 
Fig. |3| (right), the value of cone half-angle 7 toward the Sun is shown on abscissa. We observe 
no excess of events coming from the Sun. Assuming Poisson statistics for both expected 
background Nb and observed events N j, s we obtain 90% C.L. upper limits on additional 
number of signal events iV| as follows [ 32 1 : 



e _ (i V B +jv| ) y^ n °»° ( N B+^l°r 

Corresponding upper limits _/V|° for each integer value of cone half-angles in interval 1 -j- 30° 
are also shown in Fig. |3] (right). Further we estimate the values of cone half-angle 7 where 
90% of expected signal events are collected depending on the mass of dark matter particle. 
And finally, considering in details transport of high energy neutrinos from the Sun to the 
Baksan detector we translate above experimental results into the limits on muon fluxes from 
the Sun direction, dark matter annihilation rates and elastic scattering cross section of DM 
particle off proton. 

4. Transport of neutrinos produced from dark matter annihilations in the 
Sun to the Earth 

In this Section we describe numerical procedure which we use to simulate muon signal in the 
BUST resulting from dark matter annihilations in the Sun. For this purpose the well known 



WimpSim package is available 25, 26]. However, to have more flexibility we use our own C 
code for MC simulation of neutrino propagation from the center of the Sun to the level of the 
Baksan detector. Our general scheme of MC simulations is very similar to that of WimpSim 
which is very well described in Ref. [26|. They only differ in several details to be discussed in 
a moment. 

In general, the initial spectrum of (anti)neutrino from dark matter annihilations in the 
Sun is a mixture of inclusive (anti) neutrino spectra from different annihilation channels = 

dNl . dN I . 

2J^ Bi-jj^ 2 -, where is the energy spectrum of j-th type of (anti)neutrino produced in i-th 
type of annihilation channel. The quantities Bi are the branching ratios of different channels 
which can be calculated given a model for dark matter candidate. Below we simulate neutrino 
propagation for several values of masses of dark matter particle from 10 GeV to 1 TeV for 
three dominant annihilation channels: bb, t~t + or W + W~ . The channel bb represents an 
example of "soft" spectrum while W + W~ and t + t~ are examples of "hard" spectra. 

We do not calculate the neutrino spectra at the production point in the Sun for each 
annihilation channel by ourselves but take them from Ref. |33|]. We check that they coincide 
with a good accuracy with those obtained with the help of WimpSim package except for small 
deviations in spectra for v T and v T . During the propagation of neutrinos from the center of the 
Sun to the level of the detector we take into account three-flavor neutrino oscillations in the 
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vacuum, the Sun and Earth including matter effects, absorption of neutrinos due to charged 
current (CC) interactions with the interior of the Sun and Earth and the change of the energy 
spectrum due to neutrino neutral current (NC) elastic scattering and r-regeneration. 

Neutrino oscillations in 3x3 scheme are implemented according to the algorithm described 
in Refs. [j34| , |35[| . According to this scheme evolution operator analytically calculated for the 
case of constant medium density is applied to the neutrino wave function. For the case of 
varying density of the Sun one can subdivide the neutrino path into sufficiently small intervals. 
Considering the medium density inside each interval as a constant the total evolution operator 
can be approximated as product of the evolution operators corresponding to the sequence of 
intervals. We use solar model of Ref. [36|. For parameters of neutrino mixing matrix we use 
the following values 

Amli = 7.62 • 10" 5 eV 2 , |Am| x | = 2.55 • 10~ 3 eV 2 , (4.1) 
S C p = 0, sin 2 6>i2 = 0.32, sin 2 <9 23 = 0.49, sin 2 9 13 = 0.026, 

which lie within 2a range of experimentally allowed values assuming normal hierarchy for 



neutrino masses [37]. 



Our code follows straight trajectories of individual neutrinos produced at some point near 
the center of the Sun which is simulated according to the distribution p3|, 38] 



n(r) = n ex.p(-r 2 /Rl M ), where Rdm = O.Ol^VlOO GeV/m DM . (4.2) 
To propagate neutrinos from the production point we calculate the interaction length as 

Lmt{Ev) = n n (<Tgc {Ev) + a NC {Ev)) + np (aCC {Eu ) + a*C(E u )) ' (43) 

which depends on neutron and proton number densities in the Sun n n and n p , respectively, and 
the corresponding CC and NC neutrino cross sections on nucleons and . Neutrino- 
nucleon DIS total cross sections and corresponding energy and angular distributions are cal- 



culated according to the expressions presented in Ref. [39] using CTEQ6 [40] parton distri- 
bution functions. Inside the Sun, we choose the step size 5r as the smallest between Rq/F 
and L int /F where F = 300 for E v > 100 GeV and F = 3000 for E v < 100 GeV. Smaller 
steps at lower energies are required to take into consideration smaller oscillation length. We 
check that total accuracy in neutrino fluxes obtained using these step sizes is better then 1%. 
Then, we calculate the probability for neutrino to have an interaction during this small step 
Sr: P = 1 — exp(—5r /Li n t). After that we simulate the interaction process and if it takes 
place we determine its type (CC or NC, proton or neutron) with the probabilities according to 
corresponding cross sections and proton/neutron densities. When calculating CC neutrino- 
nucleon cross sections we include nonzero value of r-lepton mass while electron and muon 
are considered as massless. In NC scattering process the flavor structure of neutrino wave 
function is leaved unchanged. Neutrino energy after this interaction is simulated according 
to the corresponding differential cross section. For CC interactions the flavor of neutrino also 
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has to be simulated and in the case of u e and neutrinos we just drop them from the flow. 
For u T we simulate the energy of resulting r-lepton. Then the secondary neutrinos resulting 
from r-decay are considered following procedure described in Appendix of Ref. which is 
somewhat different from WimpSim implementation of this process 1 . 

We check results obtained by our code for neutrino differential fluxes at the level of 
the Earth (1 a.u.) with that of obtained by WimpSim and the comparison is presented 
in Figs. ||, [| and |6| for three different annihilation channels and different masses of dark 
matter particle. The results agree with each other sufficiently well. The main distinctions are 
related to different initial spectra for r-neutrinos which are also shown for t + t~ and W + W~ 
annihilation channels. 

High energy muons to be detected in neutrino telescope are produced in CC interactions 
of neutrinos in the rock just below the detector. The muons lose their energies during their 
passage from the production point to the detector. The average muon energy loss is usually 
parametrized by the following formula 

(—) = -(a(E) + (3(E)E)P, (4-4) 

where p is the density of the rock and the quantity a(E) corresponds to ionization energy loss 
while /3(E) accounts for bremstrahlung, photonuclear interactions and pair production. We 
take into account the energy dependence of a(E) and /3(E) to simulate the muon propagation 
in the rock. In the interval 1 -j- 1000 GeV of muon energies these coefficients change p2| from 
1.8 • 10~ 3 to 2.7 • 10~ 3 GeV-cm 2 /g for a(E) and from 0.8 • 10~ 6 to 3.9 • 10~ 6 cm 2 /g for (3(E). 
Having rather moderate energy dependence in a(E) and /3(E) in many calculations constant 
values for these quantities are used, for instance the values a(E) ~ 2.2 • 10 -3 GeV-cm 2 /g and 
/3(E) ~ 4.4 • 10 -6 cm 2 /g for the rock are used in WimpSim. One could expect, that for low 
energies E^ < 20 GeV neglecting energy dependence for these coefficients leads to shorter 
muon ranges and some underestimate of muon flux. At the same time at high energies as 
En > 200 -j- 300 GeV it should result in longer muon ranges and weak overestimate of muon 
flux. In Fig. ^ we compare the both resulting muon energy spectra d^^/dE^ calculated 
without and with energy dependence in ot(E) and /3(E) using values tabulated in Ref. [42]. 
We see that the obtained results confirm expectations discussed above. Generally, we find that 
the correction due to the energy dependence of a(E) and /3(E) can reach 18% to the resulting 
muon flux in the rock from bb annihilation channel for ?tidm = 10 GeV and about 2-3% to 
one from annihilation into W + W~ for th-dm = 1000 GeV. So, the effect is more pronounced 
at low energies. In Fig. || (left) we show the ratio of muon ranges R fl (E, E tn ) / R approx (E , E tn ) 
where Rp,(E, E t h) is calculated from Eq. (|4.4|) taking into account energy dependence in a(E) 
and /3(E) while Rf prox (E,E th ) is obtained for constant values a(E) » 2.2 • 10~ 3 GeV-cm 2 /g 
and /3(E) « 4.4 • 10" 6 cm 2 /g. As a cross-check we also verify our simulations of muon flux 
by comparison with an analytical expressions available for constant values of a(E) and /3(E) 
(see e.g. Ref. @). 

1 Secondary neutrinos contribute only a small part of total neutrino flux and the differences between the 
two approaches are not meaningfull. We thank Joakim Edsjo for correspondence on this point. 
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The last step of calculation is to estimate of cone half-angle which contain 90% of events 
of upgoing muons expected from dark matter annihilations in the Sun for three different 
annihilation channels bb, t + t~ and W + W~ . In this calculation we simulate the direction 
and energy of muon produced after CC interactions of neutrino in the rock, muon energy 
losses during their passage on the way to the detector and also take into account slight 
widening of muon angle distribution due to multiple Coulomb scattering [24]. The obtained 
dependence of half-cone angles on mass of dark matter particle is shown in Fig. ||| (right) and 
in Table. [I] (see next Section). 

5. Results 

The upper limits on the number of signal events iV| obtained from the results of the Baksan 
Underground Scintillator Telescope (see Section 3) can be used to put the 90% upper limit 
on the flux ^ of muons with energies higher then 1 GeV from the dark matter annihilation 
rate. This recalculation can be schematically written as 

AT 90 

where T is the live time, e = 0.9 accounts for the fact that we collect only 90% of total signal 
events inside corresponding cone half-angle and S e ft is effective area of Baksan telescope in 
the direction of the expected signal which is averaged with Sun directions and also includes 
the muon detection efficiency Q. Correspondent effective exposure towards the Sun 
amounts to 0.49 x 10 15 cm 2 s. The muon flux is related to the dark matter annihilation rate 
Fa in the Sun as follows 

F A f m DM 

= — ^ x > / dE v .P(E u .,E th ) — VJ -, (5.2) 

where R = 1 a.u., dN Vj / dE Uj is the spectrum of j-th type neutrino at the Sun surface, 
P(E Uj , E t h) is the probability to detect muon with energy E^ > E t h = 1 GeV from Uj 
neutrino with energy E Uj and the sum goes over all types of neutrino and antineutrino. The 
probability P(E Uj , E t h) has been calculated as described in the Section 4 for several masses 
of dark matter particles and annihilation channels bb, t + t~ and W~*~W~. 

In the case when the processes of dark matter capture and annihilations in the Sun reach 
exact equilibrium during time intervals much less than age of the solar system, 4.5 Giga years, 
the annihilation rate Ta should be equal to a half of capture rate determined by scattering 
cross sections of dark matter particles off solar matter. In that case, annihilation rate can be 
divided into two parts corresponding to either spin-dependent (SD) or spin-independent (SI) 
type of contribution to neutralino-nuclei interactions as follows 

Ta = Fa(<7si) + Fa(osd), (5.3) 
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where Ta(ctsi) and Fa(&sd) are the parts of equilibrium annihilation rate determined by 
either SI or SD interactions. Using the upper limits on the dark matter annihilation rate in 
the Sun and hypothesis about exact equilibrium between capture and annihilation processes 
we can set conservative upper limits on SD and SI elastic cross section of dark matter particles 
on proton p5| . To reach this goal we use recalculation procedure described in where we 
obtained corresponding coefficients in the following expressions 

4 P n PLm (m x ) = X SD (m x ) • T U /^ m (m x ) (5.4) 
a U s ? pLim (m x ) = X SI (m x ) • T U /^ m (m x ) (5.5) 

as functions of neutralino mass m x . In obtaining X SD (m x ) and X SI (m x ) we assumed 
0.3 GeV/cm 3 for the local density of dark matter and 270 km/s for root-mean-square of 
the velocity dispersion. For other details we refer reader to Ref. p6| . 

The results for upper limits on quantities discussed above are collected in Table [l] for 
several masses of dark matter and for chosen annihilation channels. Here we present obtained 
half-cone angles 7 and the upper limits at 90% C.L. on following quantities: the number of 
signal events N% , muon flux from the Sun in the rock, annihilation rate in the Sun Fa, the 
spin-dependent a~P an( \ spin-independent a^L scattering cross sections of dark matter particle 
off proton and also upper limits on neutrino fluxes <!>„ integrated from energy E u = 1 GeV. 

In Fig. H (left) are shown the 90% C.L. upper limits on muon fluxes from the Baksan's 
data for the annihilation channels W + W~, bb and t + t~ along with corresponding limits of 



SuperKamiokande [|T|, IceCube |J and ANTARES g3|, gj. T yP ica lly, for indirect searches 
with neutrino telescopes the annihilation channel t + t~ are not considered separately from 
W + W~ to set limits on dark matter properties. However, the difference between neutrino 
spectra from DM annihilation into t + t~ and W + W~ plays an important role. In spite of 
the fact that the limits on muon fluxes are approximately the same for W + W~ and t + t~, 
the limits on SD cross section (and on the corresponding annihilation rate) are quite different 
for these two channels as one can see in Figs. 10. The explanation of this fact resides partly 
in the effect of neutrino oscillations: the enhancement due to this effect for t + t~ channel is 
considerably larger than for the case of W + W~ channel. This is illustrated in Fig. || (right) 
where the ratio of muon fluxes with and without oscillation is shown for three annihilation 
channels depending on mass of dark matter. Another important effect is a breaking "democ- 
racy" in average number of high energy neutrinos per one act of dark matter annihilations: 
for t + t~ pairs it is larger in a factor of 2 4- 3 than for W + W + or bb channels. It can be seen 
in Fig. H(right) that neutrino oscillations result in a decrease of muon flux from bb annihila- 
tion channel, which is expected to be dominating at lighter neutralino masses. At the same 
time the muon flux from W + W~ annihilations is enhanced. And finally, using t + t - channel 
allows to extend the upper limits on SD cross section for "hard" type of neutrino spectra to 
lower masses of dark matter. 



The Baksan 90% C.L. upper limits on SD cross section are shown Fig. 10 in comparison 
with experimental results from direct and indirect DM searches. Here we show recent lim- 
its of SuperKamiokande 0, IceCube @, ANTARES gj|, 0> PICASSO @, KIMS @, 
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SIMPLE [§g], CMS H and ATLAS |T| collaborations. A signal range of DAMA |§ [jj 
experiment is also depicted. Strictly speaking, the indirect limits on SD and SI cross sections 
are not fully model independent. In particular, the limits of neutrino telescopes strongly rely 
on hypothesis about equilibrium between capture and annihilation processes with dark mat- 
ter particles in the Sun that may not be satisfied in some models. Impressive collider bounds 
by ATLAS [|5l]] and CMS []52| collaborations on dark matter interactions were obtained by a 
recalculation from direct limits on cross sections of dark matter pair production in association 
with jet 2 in pp-collisions. The recalculation [55, [56|, 57] is based on the fact that the diagrams 
describing the process of dark matter pair production in pp-collisions and elastic scattering 
process of DM off proton are related by crossing symmetry. However, these procedure and 
corresponding indirect bounds on SI and SD cross sections are model dependent because of 
different kinematics of these processes. The collider bounds especially weaken in the models 
when dark matter interacts with quarks via light mediators [55, [56|, [5?]] which is explained by 
completely different scaling of dark matter pair production and dark matter -nucleon scatter- 
ing cross sections with c.o.m. energy and mass of mediator. Therefore, all these indirect limits 
on SI and SD elastic cross sections of dark matter particle off proton should be interpreted 
with care. Finally, we stress that given a particular model different kinds of experiments on 
DM searches give complementary bounds on the properties of dark matter particles. 



6. Summary and conclusions 

We have performed updated analysis of indirect dark matter search with the Baksan Un- 
derground Scintillator Telescope data for 24.12 years of live time. In search for an excess 
of upward going muons in the direction toward the Sun we do not observe any significant 
deviation from expected atmospheric background. We have presented the 90% C.L. upper 
limits on spin-dependent and spin-independent elastic cross section of dark matter on proton 
assuming particular annihilation channels bb, W + W~ and t + t~. These limits have been 
derived from the limits on muon fluxes and annihilation rates for dark matter masses in the 
interval 10 4- 1000 GeV. The best value for the upper limit on SD elastic cross section is 
about 3 x 10~ 4 picobarn for masses of dark matter particle within 100 -f- 200 GeV, that is 
comparable with those presented by other operating neutrino telescopes. At the same time in 
the range of low masses of dark matter particle 10 -j- 50 GeV the Baksan limits on SD elastic 
cross section in the case of t + t~ annihilation channel are the most stringent to date among 
the results obtained by neutrino telescopes. 
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Table 1: Half-cone angles 7, upper 90% limits on the number of signal events iV| , the muon flux 
<J> M , the dark matter annihilation rate in the Sun Fa, the dark matter -proton spin-dependent cr^f 
and spin- independent scattering cross sections and neutrino fluxes $„. 
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Figure 1: Left: Measured rate of upward through going muons at the BUST since December 1978. 
Right: Distribution of measured events in angle ^-Sun between incoming muon events at the BUST 
and the Sun position. The direction of the Sun corresponds to cos^-Sun = 1. 
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Figure 2: Left: Zenith distributions of the Sun during the BUST runtime: observed (red points with 
error bars) and calculated using the Positional Astronomy Library (blue). Right: Zenith distribution 
for asymmetry between the number of neutrino events coming at time when the Sun is under 

(night, N) or above (day, D) the horizon. 




Figure 3: Left: Distributions of cos^-sun for events toward true Sun position below horizon 
(red) in comparison with measured background (green) obtained from the distributions for cos "f^-sun 
averaged over six shifted Sun positions. Right: The 90% C.L. upper limits on additional number of 
upward going muon events in the direction toward the Sun (blue) in cone half-angle 7. Measured (red) 
and background (green) events are also shown depending on the value of 7. 
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Figure 4: The energy spectra of neutrino at production point and at the distance R =1 a.u. 
from the Sun obtained with the help of WimpSim package and by our calculations for bb annihilation 
channels and mDM = 10 GeV (left) and 100 GeV (right). 




Figure 5: The energy spectra of v r neutrino at production point and neutrino at the distance 
R = 1 a.u. from the Sun obtained with the help of WimpSim package and by our calculations for 
t + t~ annihilation channels and todm — 10 GeV (left) and 100 GeV (right). 




Figure 6: The energy spectra of v r neutrino at production point and neutrino at the distance 
R = 1 a.u. from the Sun obtained with the help of WimpSim package and by our calculations for 
W+W~ annihilation channels and m DM = 100 GeV (left) and 1000 GeV (right). 
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Figure 7: The comparison between muon energy spectra obtained taking into account energy 
dependence for a(E) and (3(E) (red, solid line) and using constant values a(E) s» 2.2- 10 -3 GeV-cm 2 /g 
and f3{E) w 4.4 • 1(T 6 cm 2 /g (blue, dotted line). 




Figure 8: Left: The ratio of muon ranges calculated using full energy dependence for a(E) and P(E) 
and using constant values a(E) s» 2.2 • GeV-cm 2 /g and (3{E) w 4.4 • cm 2 /g. Right: Cone 
half-angle 7 as a function of the mass ttidm for three annihilations channels 66, t + t~ and W + W~ . 
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Figure 9: The Baksan limits on muon flux from the dark matter annihilations in the Sun in 
comparison with other experimental results (left); ratios of expected muon fluxes both from DM 
neutrino and antineutrino generated in pure three annihilation branches bb quarks, or t + t~ leptons 
or W + W~ bosons in cases with and without three flavours oscillations (right). 
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Figure 10: The Baksan limits on SD elastic cross section of dark matter particle on proton in 
comparison with other experimental results. We implement the DMTools |47j] database to plot the 
results of direct searches. 
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